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ABSTRACT 


A project  to  investigate  fluid  friction  and  cavitation  phenomena  in 
unsteady  motion  by  means  of  a special  unsteady  flow  water  tunnel  of  unconven- 
tional design  has  been  under  way  in  the  Hydrodynamics  Laboratory  of  the  Massa- 
chusetts Institute  of  Technology  under  sponsorship  of  the  0M1  since  1948.  Its 
objective  is  to  obtain  for  transient  flow  conditions  through  conduits  and  around 
immersed  bodies  the  same  kinds  of  basic  information  that  have  been  determined 
for  steady  flow  cases.  The  initial  phase  of  the  program  has  involved  the  de- 
sign and  construction  of  a pilot  model  water  tunnel  to  produce  flow  in  which  the 
basic  motion  is  accelerating  or  decelerating,  and  the  development  of  suitable 
instrumentation  for  measuring  and  recording  unsteady  velocities  and  pressures. 

The  tunnel  is  of  the  blow-down  type  and  consists  essentially  of  two 
tanks  20  in.  in  diameter  and  6 ft.  8 in.  long  mounted  vertically,  one  above 
the  other,  and  connected  by  a tube,  or  working  section,  of  one  inch  inside 
diameter.  Flow  from  the  upper  tank  passes  through  guide  vanes  and  an  entrance 
nozzle,  through  the  one-inrh  tube  which  extends  inside  the  lower  tank,  and 
through  a quick-opening  valve.  High  frequency  response  electronic  cells  were 
developed  for  measuring  differential  water  pressures  directly.  These  cells 
were  employed  in  experiments  to  determine  loss  of  head  along  the  test  section 
and  for  measuring  the  instantaneous  average  velocity. 

The  velocity  and  acceleration  of  the  flow  through  the  tube  are 
governed  by  the  difference  in  air  pressure  above  the  water  surfaces  in  the  upper 
and  lower  tanks.  A servo  control  system  has  bean  designed  and  some  components 
completed,  but  for  initial  experiments,  control  was  effected  by  hana  operation 
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of  throttle  valves  to  admit  high-pressure  air  to  the  upper  tank  and  to  maintain 
constant  pressure  in  the  lower  tank. 

The  loss  of  head  due  to  fluid  friction  in  turbulent  flow  was  deter- 
mined experimentally  for  steady  and  unsteady  flow  in  a 1M  diameter  smooth  brass 
tube  9 ft.  long.  Measurements  were  made  in  a section  where  the  boundary  layer 
was  fully  developed.  Unsteady  flow  runs  were  made  with  average  velocities  of 
15  to  72  ft/sec  and  averago  accelerations  from  zero  to  35  ft/sec2. 

Results  of  the  experiments  to  date  indicate  that  there  is  no  appre- 
ciable difference  between  the  friction  factors  referred  to  the  instantaneous 
velocity  for  steady  and  unsteady  flow  as  long  as  the  rate  of  acceleration  is 
constant  or  changes  slowly.  For  sudden  impulses  from  steady  flow  to  accelerated 
flow,  the  measured  values  of  the  wall  shearing  stress,  "C  • below  those  for 

steady  flow  at  the  same  Instantaneous  average  velocities.0  This  indicates  a 
change  in  the  distribution  of  internal  shear  stresses  and  velocity  distribution 
from  the  steady  flow  distributions  at  the  same  velocity.  These  results  for 
sudden  impulses,  although  obtained  consistently,  are  considered  tentative  pending 
accumulation  of  additional  test  data. 

The  experimental  work  performed  to  date  is  but  the  initial  phase  of 
a long-range  program  and  marks  the  initial  successful  operation  of  the  tunnel 
and  instrumentation.  The  long-range  program  is  to  construct  a much  larger 
tunnel  with  a test  section  of  sufficient  diameter  to  permit  studies  on  sub- 
merged bodies  for  cavitation  and  hydrodynamic  forces  in  unsteady  flow.  Work 
wiil  be  continued  on  the  present  tunnel,  not  only  with  the  hope  of  obtaining 
more  fundamental  knowledge  in  hydromechanics  of  unsteady  flow,  but  with  the 
aim  of  developing  an  adequate  control  system  and  determining  optimum  design 
data  for  the  larger  tunnel. 
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I INTRODUCTION 


All  known  water  tunnels  constructed  prior  to  the  one  described  herein 
have  been  for  steady  flow.  While  they  have  contributed  tramendouslv  to  the 
knowledge  of  steady  flow  fluid  friction  and  cavitation  phenomena  ( 1; , they  have 
not  been  capable  of  yielding  comparable  information  for  unsteady  flow.  Xn  1948 
the  basic  features  of  a special  "unsteady  flow  water  tunnel"  analogous  in  some 
ways  to  the  blow-down  tunnels  for  air  were  conceived  and  analysed  at  the  Hydro- 
dynamics laboratory  of  the  SJaaaashuaetta  Institute  of  Technology.  Under  the 
sponsorship  of  the  Office  of  Naval  Research,  a project  to  build  a water  tunnel 
cf  unconventional  design  which  would  be  devoted  primarily  to  studies  of  unsteady 
flow  was  initiated.  The  program  was  also  to  include  the  development  of  instru- 
ments for  measuring  Instantaneous  velocities  and  pressures. 

The  overall  objectives  of  the  program,  as  set  forth  in  an  earlier 
report  (2),  are  to  lnveatigate  and  correlate  the  characteristics  of  fluid  fric- 
tion and  cavitation  for  steady  and  unsteady  flow  through  conduits  and  past  sub- 
merged bodies.  This  is  to  Include  an  investigation  of  the  characteristics  of 
boundary  layer  flow  and  an  attmapt  to  establish  whether  the  boundary  flow  and 
its  relation  to  fluid  friction  differ  from  the  laws  established  for  steady  flow. 
The  inception  and  growth  of  cavitation  as  affected  by  local  acceleration  and 
by  presence  of  dissolved  gases  ie  to  be  investigated  also. 

The  design  of  the  tunnel  and  experiences  in  developing  Instruments  were 
described  in  an  earlier  report.  Design  considerations  led  to  the  construction  of 
« pilot  model  blowdown  type  tunnel  having  a 1"  diameter  test  section.  Tunnel 
costs  Increase  rapidly  with  sise,  but  the  instrumentation  for  a small  tunnel  can 
be  etrrieH  «y§r  to  t larger  cne  with  little  alteration.  This  report  covert  the 
first  period  of  operation  of  tha  tunnel  and  contains  an  account  of  experimental 
studies  of  fluid  frlotlon  due  to  accelerated  flow  in  a conduit.  An  olaetronic 
differential  oell  is  described  whloh  is  used  to  "instantsnoously"  measure  average 
velocities  and  pressure  drope. 

The  primary  objaotlva  of  tho  present  investigation  la  to  employ  the 
tunnel  to  determine  whether  or  not  the  pretence  of  an  acceleration  affaota  the 
wall  shear  strata  in  flow  through  • olroular  conduit  of  uniform  diameter.  It 
was  flret  proposed  to  obtain  overall  affects,  and,  If  these  war*  pronounced,  to 
then  conduct  more  detailed  experiments. 

Only  fully  developed  turbulent  flow  at  distances  of  not  lata  than  60 
diameters  from  ths  entrance  to  tha  teat  section  were  considered  in  this  study. 
Staady  flow  runs  wart  made  with  Reynolds  numbers  from  75,000  to  750,000  (svarsgt 
velocities  about  9 to  65  ft/sec) • Unsteady  flow  runs  wera  made  with  average 
velocities  from, about  15  to  75  it/ aac  and  average  accalarations  from  aaro  to 
about  35  ft/aae*1. 

A discussion  of  soma  features  of  unsteady  flow  in  tube*  is  presented  in 
Section  II  which  follows.  The  experimental  equipment  and  procedures  are  described 
in  Section  III  and  in  the  Appendix.  Tho  results  of  teat  measurements  and  dis- 
cussion of  them  are  given  in  Section  IV. 
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II  SOME  FEATURES  OF  UNSTEADY  UNIFORM  FLOW  IN  TUBES 


a.  Acceleration  In  Steady  ys  yaatiadv  Fluid  Motion* 

Fluid  motion  is  classified  as  steady  or  unsteady,  uniform  or  non- 
uniform.  Unsteadiness  means  a variation  of  the  local  velocity  with  time.  Ncn- 
uniformity  means  a variation  of  the  local  velocity  from  point  to  point  along  a 
streamline.  At  the  entrance  to  a tube  from  a reservoir,  non-uniform  flow  Is 
obtained  as  the  boundary  layer  thickness  growe  radially.  Uniform  flow  results 
after  some  distance  when  the  velocity  profile  becomes  fully  developed  by  the 
radial  convergence  of  the  boundary  layer. 

The  study  of  unsteady  motion  is  one  of  the  approaches  for  examining 
the  effects  of  acceleration.  A virtue  of  using  unsteady  motion  lies  in  the 
ability  to  cause  local  acceleration  independently  of  velocity  and  boundary 
oeometry.  Only  by  this  method  can  the  simple  case  of  acceleration  of  a uni- 
form flow  be  studied.  With  steady  motion,  acceleration  is  possible  only  with 
non-uniform  flows.  In  thsse  cases,  convective  acceleration  of  a fluid  element 
as  it  moves  through  the  fluid-filled  space  depends  upon  the  magnitude  and  di- 
rection of  the  velocity  and  the  shape  of  the  boundary. 

There  have  been  previous  studies  of  accelerated  flow  in  tubes.  First, 
the  steady  non-uniform  flow  at  the  entrance  to  a tubs  from  a reservoir,  wherein 
the  circular  core  is  accelerated  as  the  boundary  layer  grows  along  the  tube, 
has  been  subject  to  several  investigations.  These  include  the  early  study  by 
Boussinesq  (4)  in  1890  and  by  Schiller  (5)  and  later  by  Langhaar  (6),  and 
others.  Second,  cases  of  unsteady  uniform  flow  have  been  examined  including 
analytical  studies  of  oscillatory  motions  in  U-tubes  and  long  straight  tubes  (8), 
(9),  and  (10),  and  laminar  flow  in  smooth  round  tubes  (ll).  As  already  noted, 
the  unsteady  measurements  reported  here  are  for  uniform  turbulent  motion  in  a 
smooth  round  tube. 

B.  Kind  of  Flow  Under  Investigation 

As  a basis  for  examining  the  kind  of  flow  to  be  expected,  use  is  made 
of  the  boundary-layer  equation  developed  by  Prandtl  ( 12)  and  Von  Kaiman  (3) 
which  is  obtained  by  equating  the  rate  of  change  of  momentum  in  the  boundary 
layer  to  the  sum  of  the  pressure  drop  and  shear  stress  forces* 

For  simplicity,  consider  the  case  of  two-dimensional  flow  between 
parallel  flat  plates  described  by  the  diagrams  below.  The  momentum  equation, 
for  flow  with  a pressure  gradient  in  the  direction  of  motion  only,  can  be 
written 

h h h 

y p d? dy  + fx  jf  u*dy  •/,ui  f d? dy = ■ h • % 

0 0 o 

where  x is  distance  aiong  a plate,  y is  distance  normal  to  a plate,  u.  is  the 
velocity  of  the  flow  at  y = h,  u is  the  local  velocity  within  the  boundary 


layer  at  a distance  y from  one  bounding  plate,  h is  half  the  distance  between  the 
two  plates  (and  is  independent  of  time  and  distance  x) • p is  local  pressure  in- 
tensity, and  T'  is  local  boundary  shear  stress._  This  equation  is  customarily 
applied  to  turbulent  motion  by  replacing  u with  u,  the  temporal  mean  velocity. 


This  equation  indicates  a relationship  between  pressure  gradients, 
velocity  gradients,  local  acceleration,  convective  acceleration,  and  wall  shear 
stress.  It  is  to  be  noted  that  this  equation  does  not  account  for  rotational 
components  or  for  a pressure  gradient  normal  to  the  wall.  Rotational  components 
normal  to  the  wall  are  considered  negligible  in  the  experimental  setup  because 
appropriate  measures  have  been  taken  to  eliminate  such  flow  at  the  entrance  to 

the  test  tube.  In  general,  as  well  as  ^ A),  and  the  first  term  on  the  right- 

hand  side  of  Equation  [ l]  could  be  written 

/ dy 

J o 

This  means,  however,  that  if  A>»  there  is  a normal  acceleration,  i.e. , 

curvilinear  flow.  Acceleration  normal  to  the  wall  is  known  to  occur  in  this 
region  of  boundary  layer  development  near  the  entrance  to  a tube;  on  the  other 
hand,  sustained  flow  curvature  in  a long  duct  is  improbable,  and  from  conti- 
nuity considerations,  the  normal  velocity  component  becomes  a smaller  and 
smaller  percentage  of  the  longitudinal  component  as  the  boundary  layer  becomes 
fully  developed.  The  reduced  curvature  must  result  in  reduced  normal  accele- 
ration, and  hence  in  negligible  pressure  gradients  normal  to  the  wall.  This 
argument,  of  course,  assumes  no  compressibility  effects  such  as  might  lead  to 
standing  waves. 


The  work  reported  here  is  based  on  measurements  60  diameters  and 
beyond  the  entrance  to  the  tube  where  fully-developed  flow  can  be  reasonably 
expected)  hence,  dp/dy  is  taken  as  zero.  In  this  zone,  the  flow  is  uniform, 
and  at  any  instant  of  time  there  is  no  change  in  velocity  of  a filament  along  x; 
hence,  the  second  and  third  terms  of  the  left-hand  side  of  Equation  [ l]  drop  out. 
Also,  at  any  instant  of  time,  dp/dx  is  a constant  in  this  zone.  The  remaining 
term  on  the  left  side  represents  the  acceleration.  For  any  distribution  of 
acceleration  over  the  cross-section,  this  integral  is  equal  to  the  product  of 
the  average  acceleration  and  the  cross-sectional  area,  since 
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for  the  two-dimensional  flow  which  Equation  [ l]  represents. 


Equation  [i]  now  has  been  reduced  to  the  very  fundamental  Newtonian 
expression  p hi  = £"  t + Apt,  an  expression  which  might  have  been  written 
directly.  'The  development  through  application  of  the  momentum  relation,  however, 
serves  to  point  out  that  raaY  no  longer  be  the  XQ  of  steady  velocity  flow. 

C.  Wail  Shear  Stress  and  the  Distribution  of  Mean  Velocities 

In  steady  motion,  there  is  a definite  interrelation  between  the  distri- 
bution and  magnitude  of  turbulence,  the  distribution  of  the  mean  velocity,  and 
the  wall  shear  stress.  Thus  it  is  logical  to  infer  that  for  the  fully  developed 
region  under  consideration,  a modification  of  the  shape  of  the  velocity  profile 
that  might  be  caused  by  an  acceleration  would  alter  the  shear  stress  and  vice- 
versa.  Clearly,  it  is  desirable  to  study  both  the  kinematics  and  dynamics  of 
unsteady  motions  in  order  to  have  a complete  picture  of  the  acceleration  effects. 
In  the  present  study,  the  direct  measurement  of  the  velocity  distribution  was 
not  practical  because  of  the  small  size  of  the  1"  diameter  tube.  Consequently, 
it  was  possible  only  to  measure  the  change  in  pressure  gradient,  calculate  the 
wall  shear  and  infer  the  probable  change  in  velocity  distribution. 

It  is  interesting  to  consider  qualitatively  the  possible  effects  of 
acceleration  on  the  relationship  between  shear  stress  and  velocity  distribution 
with  turbulent  motion.  In  steady  uniform  flow  in  a 'conduit,  the  linearly  dis- 
tributed mean  shear  stress  is  developed  as  a consequence  of  a particular  distri- 
bution of  velocity  and  turbulence.  The  turbulence  pattern  in  the  central 
portion  of  the  stream  is  maintained  by  the  continuous  generation  and  diffusion 
of  turbulent  eddies  due  to  high  shear  forces  in  a narrow  zone  near  the  solid 
boundary.  In  fact,  this  turbulence  pattern  is  developed  by  a spreading  of  the 
boundary  influence  from  all  sides  until  a convergence  at  the  center  gives  the 
uniform  condition.  Thus,  the  establishment  of  uniform  flow  is  dependent  upon 
a time  and  space  factor.  At  conduit  transitions,  a new  regime  of  uniform  flow 
is  not  established  immediately  but  appears  after  son*  suitable  transition 
distance.  An  example  is  in  a contraction  where  the  magnitudes  of  the  velocity 
fluctuations  are  only  slightly  affected  in  passing  through  the  transition.  One 
might  conclude  that  if  an  established  uniform  motion  is  accelerated  from  one 
velocity  to  another,  some  transition  time  may  be  necessary  before  a new  es- 
tablished condition  is  effected.  In  other  words,  a time  lag  may  occur  to  allow 
the  turbulent  eddy  generating  mechanism  to  adjust  to  new  /conditions  and  to  re- 
inforce and  rearrange  the  turbulent  pattern.  With  the  existence  of  auch  a time 
lag,  the  velocity  distribution  and  hence  the  wall  shear  stress  may  not  be  the 
same  during  the  transition  ptuiod  as  for  steady  flow  at  the  same  instantaneous 
velocity.  Any  difference  would  depend  upon  the  magnitude  of  the  time  lag  and 
the  absolute  time  increment  that  such  a flow  readjustment  takes.  Experiments 
should  indicate  whether  or  not  appreciable  differences  exist. 

Next  it  is  interesting  to  consider  the  possible  effects  of  accelera- 
tion on  the  mean  local  velocity  distribution  and  the  wall  shear  that  might  be 
encountered  during  this  transition  period.  Kinematically  there  are  several 
possibilities. 

Figure  1 shows  a steady  velocity  profile  in  (a)  arxi  three  profiles  (b), 
(c),  and  (d)  representing  alternate  assumptions  concerning  the  effect  of  an  acr 
celeration  acting  for  a time  At. 
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Fig,  44  Percentage  Change  in  Friction  Factor  for  Cnange 
in  Velocity  Distribution 


In  (b)  each  filament  is  accelerated  by  the  same  amount.  Hence,  pro- 
file (b)  will  have  the  identical  shape  as  that  in  (a),  but  the  form  will  have 
changed  and  the  velocity  will  be  more  uniformly  distributed  over  the  cross 
section.  If  filaments  close  to  the  wall  are  accelerated  as  much  as  those  at 
the  center  while  particles  at  the  wall  remain  at  rest,  the  velocity  gradients 
at  the  wall  and  the  wall  shear  stress  should  be  higher  than  for  the  same 
average  velocity  with  steady  flow. 


The  profile  in  ( c)  has  the  same  form  as  that  in  ( a) , but  not  the  same 
shape.  For  example,  if  y.  exponential  equation  describes  profile  (a),  the  same 
equation  describes  pro.jJe  (c). 
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Thus,  in  this  case,  the  acceleration  is  distributed  over  the  cross  section  in  the 
same  manner  as  the  mean  local  velocity. 

As  a third  case,  let  it  bo  assumed  that  the  central  portion  is  accele- 
rated more  than  in  the  second  case,  causing  the  velocity  profile  to  become  more 
"peaked"  with  subsequent  relatively  lower  velocity  gradients  at  the  wall.  This 
case  is  Illustrated  by  Figure  1 (d)  wherein  the  profile  near  the  wall  at 
(t  + At)  is  chcv^.  approximately  roincident  with  that  at  time  (t)  while  the 
average  velocity  has  increased.  In  this  case,  neither  the  shape  nor  the  form 
of  the  velocity  distribution  curve  has  been  preserved,  and  during  acceleration 
T corresponding  to  an  instantaneous  average  velocity  would  appear  to  be  less 
than  for  steady  flow  at  the  same  average  velocity. 

It  is  not  easy  to  decide  j priori  which  of  these  three  possibilities 
will  exist.  A time  lag  could  conceivably  lead  to  either  of  the  two  extremes 
represented  by  (b)  and  (d),  whereas  if  turbulence  generation  and  diffusion  is 
rapid  enough,  case  ( c)  would  be  realized. 

D.  Estimates  of  the  Effect_Ql Velocity  Distribution  Change  on  Shear  Stress 

It  is  interesting  to  establish  a measure  of  the  deviation  of  the 
velocity  distribution-shear  stress  relation  during  transient  conditions  from 
that  for  the  equivalent  steady  flow.  As  already  noted,  this  can  be  done  by 
comparing  experimentally  determined  shear  stresses.  A simple  estimate  of  the 
magnitude  of  the  r ffect  can  be  obtained  using  the  power  law  velocity  distribution 
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If  n is  taken  as  1/7  and  C as  8.74,  the  expression  is  valid  in  the  Reynolds  number 
range  from  3000  to  100,000  and  for  higher  ranges  m is  taken  as  1/8,  1/9,  etc.  In 
Figure  3,  pipe  friction  factor  curves  calculated  for  several  values  of  m and  C are 
compared  with  the  Karman-Prandtl  logarithmic  pipe  friction  factor  curve  over  the 
range  of  Reynolds  numbers  considered  in  this  report.  Figure  2 shows  the  velocity 
profiles  represented  by  m values  from  1/6  to  1/10. 


If  it  is  now  assumed  that  the  velocity  profile  (at  least  in  the  neighbor- 
hood  of  the  wall)  retains  an  exponential  form  during  an  acceleration  and  that  only 
the  constant  and  exponent  change,  predictions  as  to  change  in  form  can  be  made 
from  observations  of  change  in  wall  shear  stress.  Figure  4 is  a plot  of  per- 
centage change  in  pipe  friction  factor  corresponding  to  deviations  in  the  form 
of  the  velocity  distribution  profile  from  that  described  by  ra  = 1/8  and  C = 9.60 
which  is  arbitrarily  taken  as  reference  base.  For  example,  it  can  be  seen  from 
Figure  4 that  at  a velocity  of  50  ft/sec.,  a gain  in  the  friction  factor  of 
4-1/256  over  the  friction  factor  corresponding  to  steady  flow  at  the  same  velocity 
would  indicate  a steepening  of  the  profile  from  the  l/8th  to  the  1/9 th  exponen- 
tial; on  the  other  hand,  a decrease  of  about  6-1/296  in  f would  indicate  the  less 
steep  profile  described  by  the  l/7th  exponential. 

It  is  not  proposed  here  that  these  exponential  laws  actually  represent 
the  velocity  distributions  in  accelerated  flow.  They  do,  however,  serve  as  a 
means  of  qualitatively  expressing  a change  in  form  of  the  velocity  distribution 
profile  corresponding  to  a change  in  shear  stress. 


Ill  EXPSII  MENTAL  EQUIPMENT  AH)  PROCEDURES 


A.  The  Unateadv  Flaw  Water  Tunnel 

The  design  considarationa  and  conatructlon  details  of  the  tunnel  are 
set  forth  in  Reference  #2.  Th«  primary  elements  consist  of  two  20"  diameter  by 
6'  -8”  long  tanks  connected  by  a 1"  diameter  test  section  as  shown  in  Figures  6 
and  7.  Eight  piezometer  taps  are  located  along  the  test  section  at  distances  of 
0,  9 , 21,  33,  45,  60,  79*5  and  99  diameters  (inches)  from  the  entrance  as  shown 
in  Figure  7. 

Driving  force  is  supplied  by  reflating  the  air  pressure  above  the  two 
water  surfaces.  An  electronic  control  system  is  being  developed  for  the  tunnel, 
but  the  work  reported  herein  was  carried  out  by  hand  operation  of  t:irottle 
valves  to  regulate  air  pressures  in  the  two  main  tanks. 

B.  Instrumentation 

For  this  investigation  of  loss  of  head  due  to  unsteady  flow  in  con- 
duits, three  variables  are  to  be  determined. 

1.  Instantaneous  head  loss  along  the  tube 

2.  Instantaneous  velocity 

3.  Instantaneous  acceleration 

In:  steady  flow  measurements,  accuracy  is  of  prime  importance  and  time 
is  not  ordinarily  considered  a major  factor.  In  unsteady  flow  measurements, 
however,  it  is  essential  that  all  variables  be  measured  at  the  same  instant) 
hence,  it  becomes  imperative  to  have  high  frequency  response  indicators  and 
recorders. 


Diaphragm  type  pressure  cells.  Incorporating  several  types  of  electronic 
measuring  principles,  have  been  designed  and  employed  in  numerous  investigations. 
The  diaphragp  cell  exhibits  high  frequency  response  and  requires  but  a small 
amount  of  fluid  to  actuate  it.  In  general,  diaphragm  type  transducers  cannot 
be  relied  upon  to  have  an  error  of  much  less  than  l£  of  full  scale,  and  this 
sacrifice  in  accuracy  is  justifiable  only  when  the  high  frequency  response 
characteristic  of  this  type  of  transducer  is  required.  When  measuring  total 
pressures,  a 1$  error  may  often  be  completely  satisfactory,  but  when  a pressure 
difference  between  two  points  is  desired,  it  is  quite  possible  that  two  such 
errors  in  opposite  directions  in  each  of  two  gages  may  be  equivalent  to  an 
error  of  many  times  the  magnitude  of  the  actual  differential  pressure.  To  avoid 
such  difficulties,  it  was  necessary  to  develop  a gage  which  would  measure  rapidly 
changing  differential  water  pressures  directly. 

A pressure  cell  was  developed  which  measures  differential  pressures 
independently  of  total  pressures.  Four  interchangeable  diaphracyns  and  three 
interchangeable  electronic  elements  permit  measurement  of  differential  pressures 
from  less  than  0.01  to  more  than  200  ft.  of  water.  Total  pressures  employed  have 
been  up  to  200  psi.  The  natural  frequency  cf  the  cells  with  air  on  both  sides 
is  from  500  cps  for  the  softest  diaphragm,  to  6200  cps  for  the  stiffest  one.  The 
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Fig.  6:  Schematic  Section  of  Pilot  Tunnel  Circuit 
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overall  frequency  reoponse  for  the  ceils  as  applied  in  this  study  (including  water 
filled  lead  lines,  piezometer  taps,  amplifiers,  and  galvanometers)  has  been  found 
to  be  at  least  160  cpc  to  300  cps  for  three  types  of  diaphragpns  used  for  measuring 
rapidly  changing  differential  pressures.  The  softest  diaphragm  was  employed  only 
for  measuring  tank  levels.  A complete  description  of  the  differential  pressure 
cell  developed  for  this  investigstlon  is  to  be  found  in  the  Appendix. 

The  average  instantaneous  velocity  is  measured  by  utilizing  the  inlet 
nozzle  as  a venturi  section  and  accounting  for  the  acceleration  head,  and  a check 
is  provided  by  recording  the  water  level  in  the  tanks  as  a function  of  time. 
Differential  pressure  cells  are  employed  in  both  cases.  In  the  former,  of  course, 
the  cell  measures  the  differential  head  at' the  venturi.  In  the  latter,  the  cell 
measures  the  difference  between  the  pressure  in  the  bottom  of  a tank  and  the 
pressure  in  the  air  section,  i.e.  the  water  level.  An  analysis  of  the  accelera- 
tion effect  on  venturi  measurements  is  given  in  D of  this  section. 

Acceleration  is  not  measured  directly  but  is  calculated  from  the  slope 
of  the  instantaneous  velocity-time  record.  It  would  be  possible  to  employ  an 
electronic  square-rooter  and  differential  analyzer  to  obtain  this  record  directly 
from  the  venturi  differential  head. 

The  differential  heads  representing  velocity  and  pressure  drop  along 
the  tube  are  recorded  on  a Hathaway-type  S-8  B Oscillograph.  Galvanometers 
used  in  this  investigation  have  a frequency  rasponse  of  730  cps.  Sample  records 
are  shown  to  a reduced  scale  in  Figure  8.  Two  traces  represent  tank  level,  three 
represent  heed  loss  along  the  test  section,  and  one  is  the  velocity  indicated 
by  the  throat  venturi.  Time  lines  appear  on  the  chart  at  intervals  of  0.01 
second  with  each  0.1  second  appearing  darker  than  the  others.  In  Figure  8,  the 
time  lines  were  lost  in  the  reduction  and  the  onetsecond  lines  were  ruled  in 
with  a pen  for  reference. 

c.  ini  Paujrimifi 

1.  Operation  of  Tunnel  for  Steady  and  Unsteady  Flow  Runs* 

In  essence,  a run  is  made  by  blowing  water  down  from  the  upper  tank  to  the 
lower  tank  through  the  1"  diameter  test  section  and  obtaining  measurements  of 
velocity,  acceleration,  and  pressure  drops  along  the  test  section.  The  methods 
of  making  these  measurements  are  described  later  in  this  section. 

To  prepare  for  a run,  water  is  forced  from  the  lower  to  the  upper  tank 
through  the  test  section  by  opening  the  foot  valve  and  supplying  air  pressure  to 
the  space  above  the  water  sufface  in  the  lower  tank.  When  the  proper  amount  of 
water  has  flowed  into  the  upper  tank,  the  foot  valve  is  closed  and  the  air  pressure 
in  tho  lower  tank  reduced  to  about  30  to  40  psi,  depending  on  the  maximum  velocity 
of  the  test  to  be  made.  This  pressure  in  the  lower  tank  is  held  to  assure  positive 
pressure  in  the  lead  lines  to  the  pressure  gages. 

Air  is  then  admitted  to  the  upper  tank  above  the  water  surface  by  a 1" 
hand  throttle  gate  valve  to  a value  about  5 or  10  psi  greater  than  that  required 
to  produce  a predetermined  steady  velocity  in  the  test  section  with  the  positive 
pressure  in  the  lower  tank  as  noted  above. 
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At  a signal  from  ona  operator,  the  foot  valve  le  opened  and  the  rater 
begins  to  flow.  For  steady  flow  runs,  pressure  In  the  lower  tank  Is  held  con- 
stant by  bleeding  out  through  a hand  throttle  valve)  pressure  In  the  upper  tank 
Is  also  maintained  at  a constant  value  by  bleeding  In  air  from  the  high  pressure 
air  receiver  through  another  hand  throttle  valve.  Once  the  flow  Is  established, 
very  little  adjustment  Is  needed  on  either  of  these  valves  to  maintain  constant 
velocity  in  the  test  section.  Both  pressures  are  noted  on  Bourdon-type  pressure 
gages.  Actually,  the  pressure  maintained  in  the  top  tank  is  not  the  air  pressure 
but  is  the  pressure  of  water  head  plus  air  pressure  Jusz  above  the  inlet  nozzle. 

Although  the  water  level  in  the  upper  tank  fails  about  three  feet  during 
a run,  data  is  only  recorded  during  about  0.6  feet  of  this  fall.  Approximately 
1.75  feet  la  utilized  to  establish  desired  test  pressures,  0.6  feet  is  for  running, 
and  the  balance  for  turning  off  gages  and  stopping  the  flow. 

As  soon  as  the  flow  is  established  and  danger  of  water  hammer  transients 
damaging  the  differential  pressure  gages  is  over,  the  gages  are  opened  to  posi- 
tion for  measuring.  When  the  running  zone  is  reached,  the  recording  oscillograph 
is  turned  on  by  remote  control  (hand)  and  the  data  is  recorded. 

When  the  water  level  has  passed  the  running  zone,  the  pressure  gages 
are  shut  off  (again  to  avoid  ’.voter  haunter  effects)  and  the  foot  valve  is  closed. 
The  water  is  then  returned  to  the  upper  tank  for  another  run  or  for  calibrations, 
which  are  described  ir.  later  paragraphs  in  this  section. 

Unsteady  flow  runs  differ  from  the  steady  flow  runs  described  above  only 
in  the  method  of  operating  the  upper  tank  inlet  air  throttle  valve.  In  these  runs, 
the  running  zone  is  approached  with  a fairly  low  steady  velocity.  In  the  running 
zone,  this  hand  throttle  valve  is  opened  rapidly  and  the  resulting  impulse  accele- 
rates the  water  through  the  test  section.  When  a record  is  desired  which  will 
begin  and  end  with  steady  velocities,  the  pressure  is  eased  off  about  half-way 
through  the  running  zone. 

A set  of  tests  usually  consists  of  gage  calibrations  alternated  with 
runs,  and  the  usual  practice  is  to  stretch  this  to  three  runs  between  each  set 
of  calibrations  after  the  first  few  cycles.  About  four  or  five  minutes  are  re- 
quired to  prepare  for  and  make  a run,  while  the  actual  duration  of  a run  is  only 
about  20  or  30  secorvis,  and  as  little  as  3 or  4 seconds  of  this  may  be  in  the 
running  zone.  Calibrations  between  runs  require  about  ten  minutes. 

2.  Operation  of  Differential  Pressure  Cells 

Each  differential  pressure  gage  measuring  head  loss  along  the  tube  is 
fitted  with  a Christmas  tree"  loop,  consisting  of  a valve  (#4)  at  the  inlet  to 
the  high-pressure  side  of  the  diaphragn,  another  at  the  inlet  to  the  low-pressure 
side  (#2),  a by-pass  valve  (#3 ),  and  a bleed  valve  (#l) . The  purpose  of  the  first 
three  valves  is  to  prevent  shock  arising  from  opening  and  closing  the  foot  vaive 
at  the  test  section  from  reaching  the  diaphragm.  The  fourth  is  for  calibrating. 
Such  shocks  arising  from  accidental  failure  to  open  the  by-pass  have  been  observed 
to  overstress  and  dish  diaphragms  and  on  one  occasion  to  completely  rupture  one. 
Figure  9 shows  a gage  with  its  valving. 
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Fig.  IOi  Typical  Calibration  Curve  for  Differential  Pressure  Cell 


i c 


Mher,  * run  ii  to  b<  made,  the  by- pa  as  valve  is  opened  and  the  two  in- 
let valves  closed*  As  soon  as  the  quick-opening  fodtr  valve  at  the  bottom  of  the 
test  section  is  opened,  the  gage  inlet  valves  are  opened  and  the  by-pass  is 
closed  in  that  order*  This  operation  occurs  before  the  running  section  is 
reached*  As  soon  as  the  water  has  passed  the  section  where  velocities  are 
measured,  the  Christmas  tree  valves  are  operated  in  reverse  order,  i.e.  the  by- 
pass is  opened  and  the  inlet  valves  are  closed*  When  all  gages  have  tnus  uwn 
positioned,  the  foot  valve  is  closed  to  end  the  run. 

3*  Calibration  of  Differential  Pressure  Cells 

The  differential  pressure  gages  are  calibrated  frequently  during  each 
test  period*  The  usual  practice  is  to  calibrate  the  gages,  make  one  blow  down 
run,  and  to  then  repeat  the  cycle  until  the  completion  of  the  test  period  which 
normally  is  three  to  four  hours  long.  Each  calibration  period  normally  requires 
about  10  minutes,  and  the  gages  are  not  rmaoved  from  their  positions  on  the  rack* 

Ref  wring  again  to  Figure  9,  valves  #1  and  #4  are  opened  and  #2  and 
#3  are  closed.  This  allows  pressure  to  be  applied  to  the  upper  side  of  the 
diaphragn  while  the  lower  side  is  ventod*  For  convenience,  a short  length  of 
tubing  is  attached  to  valve  #1,  which  terminates  in  a small'  will  it  thezzero 
position  of  the  water  column  scale  on  the  upper  tank)  consequently,  the  lower 
side  of  the  diaphragm  is  not  at  atmospheric  pressure*  The  net  pressure  on  the 
gage  is  now  the  level  of  water  in  the  upper  tank  indicated  on  the  water  column, 
plus  the  air  pressure  above  the  water  surface. 

A water  manometer  and  a mercury  manometer  are  manifolded  to  measure  the 
air  pressure  above  the  water  surface  in  the  upper  tank.  Air  is  now  fed  into  the 
upper  tank  through  a small  needle  valve.  At  scuewhat  random  increments  of 
pressure  shown  on  one  of  the  two  manometers  (depending  on  the  magnitude  of  the 
pressure),,  the  needle  valve  is  closed  and  ope  operator  gives  a signal  to  take 
readings.  At  this  time,  a short  burst  is  run  on  the  recording  oscillograph,  the 
water  column  level  .'is  read,  and  the  appropriate  reading  is  made  of  the  air  pressure 
in  the  upper  tank.  Readings  are  recorded  on  a data  sheet.  A typical  calibra- 
tion chart  record  is  shown  in  Figure  8 c. 

D.  Computational  Procedure 

Only  those  test  rims  are  considered  where  the  preceding  and  subsequent 
calibrations  have  not  shifted.  A typical  calibration  curve  is  shown  in  Figure  10. 
Such  curves  are  plotted  by  measuring  the  distance  of  the  bursts  from  the  reference 
line  in  Figure  8 c and  correlating  each  burst  with  a data  sheet  value  of  known 
differential  pressure  applied  to  the  cell.  The  calibration  curves  show  the  posi- 
tion of  a trace  on  the  oscillograph  chart  versus  differential  pressure.  By 
scaling  a trace  on  a test  run  and  consulting  the  proper  calibration  curve,  one 
obtains  plots  of  differential  pressure  versus  time. 

A constant  velocity  tost  run  is  shown  in  Figure  8 a.  Velocity  is  ob- 
tained froc  the  venturi  differential  cell  trace  17  and  is  checked  against  the 
fall  in  tank  level  given  by  traces  13  *nd  14.  Traces  15,  16,  and  18  represent 
pressure  drops  from  60  to  99  diameters  from  the  entrance  tc  the  test  section 
respectively.  Head  loss  in  feet  per  foot  for  these  steady  flow  runs  is  averaged 
from  the  last  three  traces  and  friction  factors  are  calculated  using  these  head 
lossss  and  the  velocity  given  by  the  venturi  meter  and  level  gages. 
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In  general,  computation  of  unsteady  flow  runs  follows  that  for  steady 
runs  as  far  as  obtaining  values  of  velocity  and  head  losses  from  the  traces.  How- 
ever, certain  variations  in  the  procedure  axe  necessary. 

Design  of  the  tunnel  included  a pressure  tap  at  the  throat  of  the  inlet 
nozzle  in  the  foxm  of  an  opening  0.010"  wide  completely  around  the  periphery  of  the 
test  section  (see  Figure  7).  This  tap  permits  direct  employment  of  the  nozzle  as 
a venturi  for  steady  flow.  For  unsteady  flow,  however,  the  pressure  drop  across  the 
nozzle  indicates  both  the  change  in  velocity  head  and  the  inertia  force  due  to  the 
local  acceleration  of  the  fluid  through  the  nozzle. 

The  extra  pressure  due  to  this  inertia  effect  is  small,  and  so  was  deter- 
mined approximately  by  using  a value  calculated  for  irrotatlonal  motion  through  a 
two-dimensional  nozzle  having  the  same  boundary  geometry.  This  calculated  value 

<3vt 

was  0.18  , where  v^  is  the  instantaneous  average  velocity  in  the  test  section. 

The  maximum  correction  was  of  the  order  of  2%  for  the  lowest  velocities  and  highest 
accelerations  covered  by  the  tests  and  was  less  than  1$  for  the  maximum  velocities 
at  the  highest  accelerations. 

In  computing  unsteady  flow  runs,  the  velocity-time  curve  is  first 
plotted  to  a scale  which  allows  first  differences  of  velocity  to  be  calculated 
to  yield  acceleration.  The  slope  of  the  velocity  curve  is  obtained  both  by 
numerical  differences  and  tangents.  These  accelerations  are  then  used  with  the 
above  equation  to  correct  the  velocities  for  the  inertia  effect,  and  new 
velocity  and  acceleration  curves  are  plotted  against  time. 

The  area  under  the  adjusted  velocity  curve  is  now  calculated  and  com- 
pared to  the  change  in  level  indicated  by  the  level  gages  over  several  increments 
of  timej  however,  while  this  affords  a check,  the  level  indications  are  less 
reliable  than  the  velocity  obtained  from  the  venturi.  The  acceleration  curve  is 
now  computed  more  precisely  than  for  the  first  trial  mentioned  above,  and  the  area 
under  the  curve  is  checked  against  velocity  change  for  about  1/2  second  intervals. 

The  measured  pressure  drop  or  "head  loss"  during  the  unsteady  runs  is  re- 
duced to  feet  per  foot  and  plotted  versus  the  time  scale  of  the  run.  The  head  loss 
due  to  friction  in  steady  flow  for  the  instantaneous  velocities  computed  from  the 
constant  velocity  runs  is  also  plotted  on  this  time  scale  in  feet  per  foot.  Final- 
ly, the  acceleration  head  in  feet  per  foot  is  added  to  the  steady  state  friction 
head  loss  ordinates  to  permit  comparison  of  the  pressure  drop  computed  in  this  man- 
ner with  the  measured  pressure  drop  or  unsteady  flow  "head  loss."  (Note  that  in 
Figures  12,  13,  and  14  in  Section  IV,  crosses  are  placed  at  0.1  second  intervals 
to  designate  the  computed  unsteady  head  loss.  The  measured  head  loss  curve  is 
broken  to  accommodate  these  crosses  because  the  values  fall  so  close  to  each  other). 

In  addition  to  differences  in  methods  of  computing  steady  and  unsteady 
flow  runs,  there  is  a difference  in  the  data  used.  As  noted  on  page  17,  the 
average  head  loss  indication  from  three  spans  was  used  to  obtain  the  steady  flow 
friction  factors.  When  the  calibration  shifted  for  any  one  pressure  cell  during 
a run,  the  readings  from  that  gage  were  not  included  in  the  calculations.  It 
would  be  extrsmelv  difficult  to  exactly  duplicate  unsteady  flow  runs  by  hand  opera- 
tion of  the  throttle  valves  and  thereby  obtain  data  whi-h  could  be  averaged.  For 
this  reason,  a pattern  of  unsteady  flow  runs  was  adopter  'hich  provided  a check  on 
the  reliability  of  the  recorded  differential  pressures. 

Unsteady  flow  runs  were  made  by  acceleratinq  the  flow  from  one  steady 
velocity  to  a higher  one.  The  constant  velocity  portions  of  these  runs  provided 
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an  opportunity  to  chock  the  head  loss  readings  of  each  pressure  cell  with  the 
head  loss  data  obtained  from  more  extensive  constant  velocity  tests*  Moreover, 
completely  steady  floe  runs  were  sandwiched  between  the  unsteady  flow  runs  to 
provide  an  additional  check. 

It  was  observed  that  this  steady  flow  data  obtained  during  the  unsteady 
flow  tests  gave  a random  scattering  of  points  about  the  smooth  pipe  curve  similar 
to  the  scattering  shown  in  Figure  11.  It  was  noted,  however,  that  the  combine* 
tlon  of  piezometer  taps,  lead  lines,  pressure  cell,  and  amplification  system  on 
the  span  from  79*5  to  99  diameters  from  the  entrance  showed  considerably  less 
deviation  from  the  mean  value  (the  smooth  pipe  curve)  than  was  found  on  either  of 
the  other  spans  (60  to  99  diameters  and  60  to  79*5  diameters).  For  this  reason, 
only  the  measured  head  loss  data  from  the  span  at  79*5  to  99  diameters  from  the 
entrance  was  plotted. 

E.  Eitprlaant*]  »nH  Cnanutatlonal  Precision 

Precision  of  results  from  the  experimental  work  depends  on  the  accuracy, 
both  as  to  magnitude  and  time,  with  which  differential  pressures  are  recorded  on 
the  oscillograph  chart  and  on  the  accuracy  with  which  these  re  .*ris  can  be  reduced 
to  numerical  data. 

1.  Pressure  Indlcstionst 

It  has  been  found  (see  page  45,  Appendix)  that  the  differential  pressure 
measuring  systems  respond  to  frequencies  of  at  least  160  cps  which  indicates  there 
is  no  appreciable  time  lag  between  a pressure  change  and  its  detection.  Hysteresis 
ploys  a minor  role  inasmuch  as  diaphraga  deflections  increase  continuously  during 
a run  ( except  for  small  pressure  fluctuations  about  the  mean) . The  practice  of 
making  calibrations  bsfore  and  after  each  run,  and  of  computing  only  those  runs 
where  these  two  calibrations  are  the  same,  increases  the  reliability  of  recorded 
pressure  values. 

A large  portion  of  the  calibration  shifting  is  due  to  the  electronic 
equipment.  It  has  been  found  that  the  output  voltage  from  the  oscillator  which 
energizes  the  pressure  cells  ( see  Appendix)  is  subject  to  fluctuations.  These 
fluctuations  are  primarily  due  to  variations  in  line  voltage  and  to  vibrations 
caused  when  heavy  trucks  pass  the  laboratory.  Although  he  absdlute  variation 
is  small,  the  amplification  of  the  signal  from  the  p~o*  ve  cells  is  sufficient 
to  make  the  shifts  of  zero  position  on  the  oscillograph  chart  noticeable.  The 
practice  of  using  only  the  runs  where  th<>  calibrations  remain  constant  minimizes 
this  error. 

2.  Conversion  of  Data  from  Oscillograph  Chartsi 

Diaphragns  and  transformers  are  so  selected  and  amplification  of  the 
output  signals  so  regulated,  that  the  ranges  of  differential  pressures  covered 
by  a pressure  cell  during  a test  run  will  extend  over  approximately  20  cm.  on  the 
recording  chart.  Errors  arise  in  scaling  the  calibration  charts  shown  in  Figure 
8 c and  in  reducing  the  deflections  on  the  test  charts  (Figures  8 a and  8 b)  to 
differential  pressures. 

(a)  r?rors  in  Obtaining  Data  from  Calibration  Charts*  Three  observers 
have  found  they  read  the  distance  of  the  same  calibration  marks  from  the 
chart  refererce  line  within  0.3  mm.  Smooth  calibration  curves  are  drawn 
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through  tho  too  sots  of  calibration  points  for  aach  run  as  shown  in  Figure 
10  and  arrors  in  raadlng  chart  valuas  art*  tharaby  rsducad. 

(b)  Errors  in  Obtaining  Data  from  last  Run  Chart  si  Dua  to  tha  largs 
amplitude  of  dlffarantlal  prassura  fluctuations  about  a naan  valua  as  indi- 
catad  in  Figuras  8 a and  8 b,  tha  thraa  obsarvars  raad  tha  naan  valua  with- 
in 1.0  an  a v full  seals  and  within  .5  on  at  1/4  to  1/3  of  full  seals  deflec- 
tioA*  Although  tha  absoluta  valua  of  tha  fluctuations,  which  nay  ba  dua  to 
turbu lanes,  incraasas  with  valocity,  tha  anpli float ion  of  a signal  to  give 
full  chart  daflaction  for  diffarant  ranges  of  dlffarantlal  pressures  nakes 
tha  fluctuations  more  or  less  a function  of  chart  daflaction* 

Tha  overall  error  in  converting  data  fron  tha  oscillograph  charts  to 
numerical  data  nay  ba  as  high  as  1.3  nm  in  20  cm.  at  full  chart  daflaction  and 
0.8  mm  at  mall  deflections.  At  full  scale  this  is  equivalent  to  0.6£,  but  at 
low  prassura  readings  tha  error  may  ba  % or  more. 


IV  EXP ffil MENTAL  PROGRAM 


a.  Scare  8 X Sxpfiilflmta 

Fluid  friction  dua  to  both  steady  and  unsteady  flow  in  » l"  dianatar 
conduit  was  determined  for  turbulent  flow  in  tha  region  of  a fully  developed 
boundary  layer.  Steady  flow  runs  ware  made  with  Reynolds  numbers  from  75,000  to 
750,000  (average  velocities  about  9 to  85  ft/sac).  Unsteady  flow  runs  wore  made 
with  average  velocities  from  about  15  to  72  ft/sec  and  average  accelerations  from 
zero  to  about  35  ft/sec2 . 

As  was  explained  on  page  18  for  the  steady  flow  tests,  differential  pres- 
sures were  averaged  from  the  readings  between  taps  at  60N  and  79.5",  and  79.5”  and 
99",  and  60"  and  99"  from  the  entrance  to  the  test  section,  and  whan  the  calibra- 
tion shifted  for  any  on*  gage  (i.e.  span)  during  a run,  the  readings  from  that 
gage  ware  not  included.  On  tha  other  hand,  for  unsteady  flow  runs,  only  the  dif- 
ferential pressures  on  tha  span  from  79.5"  to  99"  from  the  entrance  to  tho  *est 
section  ware  used  because  this  data  was  the  most  consistent. 

It  should  be  emphasized  at  this  point,  however,  that  the  data  from  all 
three  spans  indicated  the  same  trends  in  all  portions  of  the  unsteady  flow  runs. 

In  runs  where  the  calibrations  of  the  gages  on  all  three  spans  were  constant,  the 
average  values  could  have  been  used,  and  the  results  would  have  been  very  close 
to  those  obtained  from  the  one  gage  alone.  A procedure  of  using  only  those  un- 
steady funs  where  all  gage  calibrations  remained  constant,  however,  would  have 
materially  reduced  the  data  available  fur  calculation. 

The  practice  of  accelerating  from  one  steady  velocity  to  another  during 
the  unsteady  flow  runs  (as  described  on  page  19)  produced  similar  patterns  of 
average  velocity  and  average  acceleration  for  all  runs.  In  order  to  obtain  the 
best  accuracy  in  scaling  data  from  the  records,  however,  amplification  of  the 
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signals  from  tha  dlffarantial  pressure  cells  was  regulated  to  give  full-6cale  chart 
deflection  for  different  portions,  or  phases  of  the  runs.  These  phases  h»ve  arbi- 
trarily been  called  the  initial  impulse  or  transition  from  steady  to  accelerated 
flow,  established  acceleration,  and  the  gradual  transition  from  accelerated  flow 
to  steady  flow. 

1.  Initial  Impulse  or  Transition  from  Steady  to  Accelerated  Flows 

This  phase  of  the  accelerated  flow  pattern  is  best  shown  in  Figure  12. 

It  includes  the  range  froo  an  established  steady  flow  it  t - 0 to  the  point  where 
the  acceleration  has  reached  its  maximum  value.  In  Figures  12  a and  12  b,  this 
is  at  1.0  seconds  and  1.1  seconds  respectively.  For  these  runs  the  ranges  of  the 
differential  pressure  cells  measuring  pressure  drop  along  the  test  section  was 
set  for  zero  to  6 ft/ft,  which  represented  the  full  maximum  range  of  linear  output 
froo  the  amplifiers  and  power  supplies.  The  test  runs  shown  in  Figures  13  and  14 
also  include  this  initial  impulse  phase,  but  the  pressure  cells  were  not  set  to  in- 
clude the  head  loss  and  velocity  measurements  associated  with  it.  In  Figure  12  a, 
the  maxlmun  acceleration  was  36  ft/sec*  and  the  maximum  rate  of  change  of  accele- 
ration was^77  ft/sec2 /sec.  In  Figure  12  b,  these  values  are  30  ft/sec2  and 
100  ft/sec2/sec  respectively. 

2.  Established  Acceleration  Phase* 

This  phase  represents  the  range  following  the  peak  acceleration  where 
acceleration  is  decreasing  or  is  nearly  constant.  Figures  12.  13,  and  14  all 
show  this  phase  to  some  degree,  but  for  the  runs  shown  in  Figure  13>  the  ranges  of 
the  pressure  cells  were  set  for  zero  to  12  ft/ft  to  cover  it  more  completely  than 
for  the  others.  The  rate  of  decrease  in  acceleration  for  this  portion  of  the  runs 
in  Figure  13  a to  13  f is  from  -2,0  ft/sec2 /sec  to  -14.4  ft/sec2 /sec.  Acceleration 
rates  were  from  zero  to  30  ft/sec2. 

3.  Gradual  Transition  from  Accelerated  to  Steady  Flow* 

This  phase  is  illustrated  in  the  runs  shown  in  Figures  14  a and  14  b. 

The  primary  value  of  these  runs  is  to  demonstrate  that  the  measured  pressure  drop 
in  the  unsteady  flow  tests  becomes  the  steady  flow  value  when  the  transition 
occurs  and  to  thereby  serve  as  a check  on  the  accuracy  of  measurements. 

B.  Presentation  of  Results 

Results  of  the  experimental  work  are  encompassed  in  Figures  11,  12,  13, 
and  14.  Figure  11  shows  the  experimental  points  determined  from  steady  flow  opera- 
tion of  the  tunnel  plotted  on  the  familiar  smooth  pipe  friction  factor  curve  and 
summarizes  the  steady  flow  experiments.  Results  of  unsteady  flow  runs  are  given 
in  Figures  12,  13,  and  14.  In  these  figures  for  unsteady  flow,  curves  are  drawn 
for  the  instantaneous  values  of  average  velocity  (u)  or  (Q/A),  average  acceleration 
(a)  over  the  cross  section,  measured  head  drop,  and  the  equivalent  steady  flow  head 
drop.  The  latter  is  calculated  using  the  instantaneous  velocities  in  the  steady 
flow  head  loss  equation.  The  sum  of  the  equivalent  steady  flow  head  loss  and 
the  acceleration  head  are  shown  by  the  small  crosses  at  0.1  second  intervals,  and 
the  curva  of  measured  instantaneous  head  drop  is  broken  to  accommodate  these 
crosses  where  the  values  nearly  coincide. 


c.  Anaiaaii  al  Bjnuiii 

1*  Steady  Flow  Runsi 


In  developing  tha  tunnel  and  instrumentation,  the  aim  waa  to  ewlve  sr. 
apparatus  which  could  readily  produce  uneteady  flow  patterns,  and  instruments  with 
high  frequency  reaponae  to  "Instantaneously”  record  the  variables.  They  should 
not  be  expected  to  compete  with  systmns  employing  constant  head  tanka  and  micro- 
manometers  for  determining  steady  flow  friction  factors.  Figure  11  shows,  however, 
that  the  steady  flow  friction  factors  form  a good  cluster  about  the  logarithmic 
smooth  pipe  friction  curve.  In  this  figure,  the  maximum  deviation  is  3*£>  but 
asst  points  fall  within  1$  to  1-1/2$  of  the  values  on  the  curve. 

2.  Unsteady  Flow  Runsi 

In  Figures  12,  13,  and  14,  the  small  crosses  at  0.1  second  intervals 
provide  a ready  means  of  comparing  the  measured  pressure  drop  along  the  tube 
with  the  values  obtained  by  combining  the  acceleration  head  and  the  steady  flow 
head  loss  based  on  the  instantaneous  average  velocity.  For  comparing  wall  shear 
stress  ( T)  values,  it  is  more  conven:  ent  to  subtract  the  acceleration  head  from 
the  measured  head  drop  curve  and  to  compare  the  difference  with  the  steady  state 
friction  head  loss  curve-.  In  the  discussions  which  follow,  the  variation  of  X 
for  unsteady  flow  from  that  for  steady  flow  was  calculated  in  this  manner. 

As  was  noted  on  page  21,  these  runs  have  been  divided  arbitrarily  into 
three  phasest  ( l)  initial  impulse,  (2)  established  acceleration,  and  (3)  gradual 
transition  from  accelerated  to  sr.eady  flow? 

(a)  Flow  with  Established  Acceleration!  The  maximum  variations  of  the 
unsteady  flow  X from  the  steady  flow  value  for  the  runs  shown  in  Figures 
12,  13,  and  14  axe  as  followsi 


Maximum  variation  of  7" 

o 


Fig. 
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Velocity 
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57 
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53 

17 

+ 

.87 

e 

39-l6c 

59 

13 

- 

.48 

f 

39-20A 

50 

13 

- 

.77 

14a 

38-6A 

45 

19 

- 

1.53 

b 

39-20B 

43 

23 

+ 

1.75 

This  tabulation  shows  a random  variation  of  the  X values  for  unsteady 
flow  from  those  for  steady  flow  which  is  within  the  variation  of  experimentally 
detexrsined  steady  flow  friction  factors  from  the  average  ( smooth  pipe  friction 
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Fig.  lli  Experimental  Steady  Flow  Friction  Factors 
Plotted  on  Smooth  Pipe  Curve 


Fig.  12a»  Unsteady  Flow  Run  - Initial  Impulse  from  Steady 
to  Accelerated  Flow 
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Fig.  12bi  Unsteady  Flow  Run  - Initial  Impulse  from  Steady 
to  Accelerated  Flow 
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Fig.  i3ai  Unsteady  Flow  Run  - Established  Acceleration  Phase 


Fig.  13bi  Unsteady  Flow  Run  - Established  Acceleration  Phase 
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Fig.  13cj  Unsteadv  Flow  Run  - Established  Acceleration  Phase 
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Fig.  13d»  Unsteady  Flow  Run  - Established  Acceleration  Phase 
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Fig.  13ft  Unsteady  Flow  Run  - Established  Acceleration  Phase 
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Fig.  14at  Unsteady  Flow  Run  - Transition  from  Accelerated  to  Steady  Flow 
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Fig.  14b»  Unsteady  Flow  Run  - Transition  from  Accelerated  to  Steady  Flow 
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curve)  in  Figure  11,  In  other  words,  the  experimental  work  indicates  that 
when  the  acceleration  la  established  and  xa  not  changing  rapidly  with  time, 
there  is  no  appreciable  difference  between  the  values  of  71  foT  unsteady 
flow  and  steady  flow.  As  was  noted  on  page  21,  the  experiments  included 
accelerations  of  from  zero  to  30  ft/sec*  and  rates  of  change  of  accelera* 
tion  from  - 2 ft/seeVsec  to  *14.4  ft/secVsec.  As  discussed  on  page  7, 
this  Indicates  that  the  fox*  of  the  acceleration  distribution  profile  must 
be  very  similar  to  the  velocity  distribution  profile. 

(b)  Initial  Impulse  Phases  The  test  runs  shown  in  Figures  13  and  14 
indicated  that  there  was  a consistent  difference  in  the  unsteady  flow  71 
from  the  steady  flow  value  during  the  initial  impulse.  It  did  not  appear 
to  be  fair  to  jud;ie  wholly  from  these  runs,  however,  because  this  difference 
occurred  at  the  lewer  portion  of  the  range  settings  of  the  differential  pres-* 
sure  cells.  As  discussed  in  Section  III  E on  page  20,  this  is  the  region  of 
maximum  error  in  converting  data  from  the  charts.  To  further  investigate 
this  region  additional  tests  were  perfozued  in  which  the  ranges  of  the  dif* 
ferential  pressure  cells  were  set  so  that  the  differential  pressures  at  the 
initial  impulse  were  approximately  at  the  middle  of  the  ranges.  Figures  12  a 
and  12  b show  the  difference  between  the  measured  unsteady  T and  the  steady 
7^  to  be  as  follows! 
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13.0 

- 5.1 

.90 

30.5 

20.7 

+ 1.3 

the  maximum  deviation  occurs  during  the  period  of  maximum  rate  of  change  of 
acceleration  rather  than  at  the  maximum  acceleration.  In  both  cases,  at  the 
maximum  acceleration  the  measured  head  drop  agrees  well  with  the  value  compu- 
ted from  the  sum  of  the  acceleration  head  and  the  steady  state  head  loss.  In 
Figure  12  a the  maximum  rate  of  change  of  acceleration  is  77  ft/sec2/sec,  and 
in  12  b it  is  100  ft/sec2/sec. 

In  turbulent  flow  in  a conduit,  the  shear  stresses  within  the  central 
core  are  primarily  due  to  turbulent  velocity  fluctuations.  As  discussed  on 
page  7,  it  is  indicated  that  a finite  time  may  be  required  for  the  velocity 
fluctuations  to  rearrange  themselves.  At  the  instant  the  accelerating  force 
is  applied,  the  motion  of  all  filaments  is  resisted  by  shear  stresses  which 
are  the  same  as  those  of  steady  flow.  As  the  velocity  increases  due  to  the 
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accelerating  force,  the  filaments  near  the  wall  are  resisted  by  Increased 
viscous  shear.  The  shear  on  the  filaments  within  the  central  core,  however, 
cannot  change  until  the  velocity  fluctuations  change.  If  the  latter  require 
a finite  time  to  change,  the  central  portion  of  the  flow  will  be  accelerated 
more  than  that  at  the  boundaries.  The  net  effect  Is  to  produce  a velocity 
gradient  at  the  wall  during  a rapid  impulse  which  is  less  than  would  be  ex- 
pected for  the  instantaneous  average  velocity  of  flow  through  the  tube.  In 
this  sense,  the  instantaneous  velocity  profile  form  can  be  said  to  become 
"more  peaked,"  and  correspond  to  the  third  possibility  (d)  discussed  on  page 
7.  As  discussed  on  pages  7 and  8,  the  exponential  forms  of  velocity  distri- 
bution can  be  used  to  illustrate  this  phenomenon.  In  Figures  12  a and  12  b, 
the  value  of  T for  unsteady  flow  during  a sudden  impulse  was  found  to  be 
up  to  Tf>  less  tnan  that  for  steady  flow  at  the  same  instantaneous  velocities. 
This  deviation  occurred  at  velocities  of  25  to  30  ft/sec,  and  Figures  3 and 
4 on  page  6 indicate  that  this  change  would  be  approximately  equivalent  to 
a change  from  m = 1/8  to  m = 1/6  in  the  exponential  velocity  distribution 
law. 


If  the  observation  that  these  deviations  occur  only  during  high  rates 
of  acceleration  is  typical,  it  would  indicate  a possible  lag  in  the  setting 
up  of  the  eddy  generating  mechanism.  Once  the  generating  mechanism  has 
"overcame  the  shock,"  turbulence  is  generated  at  a rate  sufficiently  in 
balance  with  other  conditions  to  give  no  appreciable  differences  from  steady 
motion. 

( c)  Gradual  Transition  from  Accelerated  to  Steady  Flow  The  transi- 
tions shown  in  Figures  14  a and  14  b are  gradual,  and  the  maximum  variation 
of  the  unsteady  flow  71  from  the  steady  flow  va^ue  is  0.6(1,  which  is  of 
the  same  order  of  magnitude  as  experimental  inaccuracy.  As  was  pointed  out 
on  page  21,  the  primary  purpose  of  *hese  runs  was  to  make  certain  the  pres- 
sure cells  indicated  reliable  preso-re  drops  aftsr  the  sharp  rises  in  pres- 
sures to  which  they  ware  subjected  during  the  runs. 

These  tests  were  in  no  way  counterparts  of  the  initial  impulse  phase. 
Such  a counterpart  would  entail  high  (negative)  rates  of  change  of  accele- 
ration for  periods  of  time  comparable  to  those  in  the  initial  Impulse  phase. 
It  could  also  be  produced  by  rapidly  decelerating  e steady  velocity  flow. 
(Note  that  thle  could  be  done  on  the  present  tunnel  by  arranging  to  rapidly 
introduce  high-pressure  air  to  the  lower  tank.) 

It  can  only  be  concluded  from  these  tests  that  they  further  verify  the 
results  of  the  established  acceleration  phase  runs)  namely,  that  there  is 
no  appreciable  difference  between  the  steady  flow  7^  and  the  unsteady  flow 
X vsluos  for  nearly  constant  rates  of  acceleration;  On  the  other  hand, 
if°a  steady  flow  were  decelerated  with  a high  rate  of  change  of  deceleration, 
it  is  probable  that  the  shear  stresses  in  the  central  core  due  to  velocity 
fluctuations  would  not  change  immediately.  In  this  case,  however,  the  shear 
stresses  on  each  filament  aid  the  decelerating  force,  and  the  layers  near 
the  wall  are  probably  retarded  more  rapidly  than  those  farther  from  the  wall. 
Such  action  would  lead  to  higher  velocity  gradients  at  the  wall  and  conse- 
quently higher  values  of  XQ  than  for  steady  flow  at  given  instantaneous 
velocities. 
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It  should  be  noted  for  these  observations  on  the  effect  of  rate  of 
acceleration  that  the  deviations  from  steady  flow  friction,  although  consistent, 
aib  small  percentagewise.  Furthermore,  it  is  possible  that  these  deviations 
could  be  explained  by  consistent  measuring  errors.  Thus,  while  the  instrumentation 
and  the  experimental  mni*urament  and  data  reduction  techniques  appear  reliable,  as 
was  discussed  in  Section  III  and  is  considered  further  in  the  Appendix,  these  items 
are  being  scrutinized  carefully  in  anticipation  of  additional  confirming  experi- 
ments. The  above  conclusions, consequently,  are  considered  tentative. 


V SIAMARY  OF  CONCLUSIONS 


These  experiments  with  fully  developed  turbulent  flow  in  a conduit 
indicate  the  following  variations  from  steady  flow  phenomena  based  on  the  arbi- 
trarily selected  phases  of  unsteady  flow  studied. 

1.  Established  Acceleration  Phasei 

( a)  No  appreciable  change  in  value. 

(b)  Approximately  the  same  velocity  distribution  law. 

2.  Initial  Impulse  Phase; 

The  variations  in  the  following  properties  appear  to  be  a function  of 
the  rate  of  increase  of  the  acceleration  and  not  of  the  actual  magnitude  of 
acceleration.  The  variations  from  steady  flow  noted  below  are  for  periods  of 
rapidly  increasing  rate  of  acceleration. 

( a)  X is  less  than  for  steady  flow  based  on  the  instantaneous  average 
velocity. 

(b)  Velocity  distribution  is  apparently  "more  peaked"  than  in  steady  flow. 

( c)  Shear  stress  distribution  is  not  linear.  While  shear  is  still  maxi- 
mum at  the  wall  and  zero  at  the  center,  at  intermediate  points  it 
probably  is  less  than  is  indicated  by  a linear  distribution. 

3.  Rapid  Increase  in  Rate  of  Decelerations 

This  phase  was  not  studied  experimentally,  but  it  appears  reasonable  to 
draw  the  following  conclusions  on  variations  from  steady  flow  based  on  the  results 
of  the  initial  impulse  phase  and  analytical  consideration. 

( a)  % would  be  higher  than  for  steady  flow  based  on  the  instantaneous 

average  velocity. 

(b)  Velocity  distribution  would  appear  to  have  higher  gradients  at  the 
wall  than  in  steady  flow. 

( c)  Shear  stress  would  be  maximum  at  wall  and  zero  at  center  but  probably 
higher  at  intermediate  points  than  in  a linear  distribution. 

The  observations  concerning  the  rate  of  acceleration  or  deceleration  are 
based  on  measured  differences  that  are  small  percentagewise  and  are  considered 
tentative  pending  additional  confirming  experiments. 
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APPBDIX 


DIAPHRAGM  DIFFERENTIAL  TRANSFORMER  TYPE  DIFFERENTIAL  PRESSURE  CH-LS 


A.  General 

The  diaphragm  type  differential  pressure  gage  developed  on  this  pro- 
ject incorporates  a small  differential  transformer,  which  has  light  weight 
moving  parts,  does  not  require  close  clearances,  and  can  be  operated  for  at  least 
several  months  under  water.  Its  development  is  the  outgrowth  of  studies  and  ex- 
periences reported  in  Reference  (2). 

The  operating  principle  of  the  differential  transformer  sensing  element 
can  best  be  explained  by  referring  to  Figure  15* 

The  transformer  consists  of  a small  iron  core  and  three  coils  wound 
side  by  side  on  tr.a  same  axis.  The  center  coil  Is  the  primary  winding,  and  the 
other  two  are  secondaries  connected  in  phase  opposition.  When  the  core  is  centered, 
the  same  voltage  is  induced  In  each  secondary  and  the  net  output  is  zero.  The 
core  is  then  said  to  be  at  its  null  point.  Movement  of  the  core  in  either  direc- 
tion causes  an  increase  in  the  coupling  between  the  primary  and  one  secondary  coil 
and  a decrease  between  the  primary  and  the  other  secondary  coil,  resulting  in  a 
net  voltage  output  which  is  linearly  proportional  to  the  travel  of  the  core. 

The  core  is  connected  by  a brass  rod  to  a diaphragm  as  shown  in  Figures 
16,  17,  and  18,  and  oovamentA  of  the  diaphragm  are  thus  detected  and  amplified. 

If  the  diaphragm  deflection  is  linear  with  applied  differential  pressure,  the  out- 
put signal  from  the  transformer  is  also  linear  with  this  differential  pressure. 

b.  ConstoucUofl  Pstaus . 

Figure  18  shows  the  assembled  differential  pressure  cell  in  cross  sec- 
tion. Details  of  mounting  the  diaphragp  and  field  coils  are  shown  in  Figures  16 
and  17.  Pressure  is  admitted  to  both  sides  of  the  corrugated  diaphragm  through 
ports  A and  B.  The  gage  is  normally  employed  with  the  greater  pressure  at  port  B, 
and  the  following  descriptions  are  for  this  arrangement)  however,  the  cell  perform 
satisfactorily  with  the  greater  pressure  at  port  A.  As  the  diophragpn  deflects, 
the  core  moves  longitudinally  along  the  axis  of  the  coils  and  a signal  is  genera- 
ted as  described  in  the  preceding  section. 

All  metal  parts  of  the  pressure  cell  are  made  of  brass  with  the  exception 
of  the  diaphragms,  which  are  corrugated  beryllium  copper  or  Inconel.  The  case  was 
made  heavy  to  guard  against  distortion  at  pressures  up  to  200  psi  where  the  gages 
are  sometimes  employed. 

The  corrugated  diaphragms  used  as  the  sensing  element  were  originally 
made  for  capsules  for  aeronautic  instruments  and  controls,  and  they  were  received 
with  a skirt  around  the  periphery.  This  skirt  was  utilized  in  mounting  in  the 
following  mannert  A heavy  brass  ring  was  machined  with  a groove  and  a lip,  as 
shown  in  Figure  17.  The  lip  was  machined  to  match  the  contour  of  the  diaphragm 
at  the  edge  and  so  that  the  skirt  on  the  diaphragm  had  to  be  pressed  lightly  over 
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Fig.  Ip.  Schematic  Diagram  of  Electronic  Principles  Employed 
in  Differential  Pressure  Cell 


Fig.  16s  Diaphragm,  Cage,  and  Transformer  Assembly 
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he  lip,  causing  a slight  radial  tension*  The  groove  was  then  filled  with  pieces 
of  solder,  the  diaphragm  was  held  in  pi  ace  with  a weight,  and  the  assmnbly  was 
rotated  on  a small  rotary  table  while  a small  torch  played  on  ihe  edge  of  the  ring 
to  melt  the  solder.  This  Job  could  also  be  done  in  a controlled  electric  oven. 

Soft  solder  was  used  in  the  groove  for  both  the  beryllium  copper  and  the  Inconel 
diaphragms.  In  early  models,  the  Inconel  diaphragms  were  silver  soldered,  but  it 
was  found  that  soft  solder  was  just  as  effective  as  the  silver  solder  and  was  much 
easier  to  apply. 

The  transformer  coil  housing  is  Inserted  in  the  adjusting  nut  as  shown. 

The  coil  asscnbly  fits  loosely  in  the  adjusting  nut  and  is  held  in  place  by  a 
lucite  sleeve,  which  is  a snug  fit.  This  design  permits  the  coil  assembly  to  be 
removed  without  damage  for  inspection  or  recosting  with  waterproofing  material. 

Standard  recommended  practice  is  to  support  the  coils  in  plastic  and  to 
keep  all  metal  as  far  as  possible  from  them.  In  early  stages  of  development,  hard 
rubber  was  used,  but  its  low  modulus  of  elasticity  and  high  Poisson  Ratio  caused 
it  to  distort  under  pressure.  Laboratory  experiments  disclosed  that  if  the  coil 
assembly  is  mounted  in  a brass  cup,  in  such  a way  that  the  metal  surrounding  the 
coil  assembly  is  approximately  symmetrical,  the  characteristics  of  the  transfozmer 
are  as  good  as  when  the  coil  is  surrounded  by  plastic.  The  location  of  the  colls 
is  not  critical  as  long  as  the  metal  around  them  is  approximately  symmetrical) 
however,  if  the  coil  assembly  is  raised  until  it  is  flush  with  the  top  of  the  ad- 
justing nut,  the  output  almost  vanishes.  When  the  housing  is  half-way  out  of  the 
nut,  the  output  characteristics  are  fair,  but  the  location  becomes  extremely 
critical,  and  change  in  location  of  a few  hundredths  of  an  inch  causes  the  signal 
to  vanish  entirely. 

The  adjusting  nut  and  coll  assembly  are  coupled  to  the  diaphragm  support 
ring  by  a cage.  This  arrangement  makes  all  relative  motion  between  diaphragm  and 
coils  independent  of  the  total  hydrostatic  pressure  in  the  shell  nf  the  pressure 
cell.  The  adjusting  nut  screws  in  or  out  of  the  cage  (40  threads  per  inch)  to 
permit  setting  the  transformer  for  the  desire;!  operating  range  as  described  in  a 
later  paragraph. 

Lead  wires  from  the  transformer  coils  were  coated  with  rubber  cement  and 
drawn  through  small  holes  drilled  in  a #1  rubber  laboratory  stopper  and  were  soldered 
to  an  amphenol  connector,  as  shown  in  Figure  18.  The  rubber  stopper  was  coated 
with  rubber  cement  and  fitted  into  a brass  bushing  sealed  with  "0"  rings  and  held 
in  place  by  a set  screw.  When  it  is  desired  to  change  a transformer  coil,  the  en- 
tire assembly  of  adjusting  nut,  coil,  bushing,  and  amphenol  connector  are  removed 
as  a unit.  All  such  units  axe  interchangeable  from  gage  to  gage. 

C.  Sensinn  and  Detecting  Elements 

1.  Diaphragms* 

Beryllium  copper  and  Inconel  diaphragms  exhibit  low  hysteresis  character- 
istics and  appear  to  be  the  two  best  materials  known  for  this  service.  Corrugated 
diaphragms  are  more  linear  and  do  not  develop  as  much  radial  tension  for  a given 
deflection  as  flat  diaphragms.  Two  types  of  diaphragms  supplied  by  Ranco,  Inc., of 
Columbus,  Ohio.,  were  used  in  these  gages.  The  first  is  a 15/16"  diameter  by  0.0075" 
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thick  Ir.conel  diaphragm  with  three  corrugations,  and  the  second  is  a 1-7/32"  by 

0.0043"  beryllium  copper  diaphragm  with  three  corrugations.  Each  of  these  dia- 
phragms was  made  to  serve  two  pressure  range?;  One  pressure  range  was  obtained 
by  using  their  full  diameter  as  shown  in  Figure  17.  The  second  pressure  range 
from  each  was  obtained  by  reducing  the  diameter,  it,  to  the  next  corrugation  and 
soldering  the  diaphragm  to  the  ring  at  this  diameter.  The  reduced  sizes  are 
designated  "half-siz*  Inconel,"  and  half-size  beryllium  copper"  diaphra^ns. 

2.  Differential.  Transformer  Elements! 

These  eluents  were  supplied  by  Schaevitz  Engineering  Corporation,  of 
Camden,  New  Jersey.  Three  sizes  are  used,  0.005L,  0.01CL,  and  0.020L.  The  number 
indicates  the  inches  of  travel  of  the  core  for  which  the  output  is  linear,  i.e. 
0.003"  etc.  When  ths  travel  is  greater  than  the  rated  value,  the  output  is  no 
longer  linear.  The  coil  assemblies  are  approximately  3/16"  diameter  by  3/8"  long 
with  a 1/8"  diameter  hole  for  the  core,  and  the  core  is  approximately  1/16"  dia- 
meter by  3/32"  long  with  a 0.030"  hole  for  mounting.  The  Schaevitz  Engineering 
Company  cooperated  with  this  project  in  developing  suitable  waterproofing  tech- 
niques and  now  supply  these  transformers  with  a heat-treated  coating  of  Glyptal 
Varnish  on  request. 

d.  Electrical  QaasapaaLa 

The  electrical  components  associated  with  Type  D C T Differential  Pres- 
sure Cells  consist  of  the  following  elements! 

1.  Oscillator! 

A Hewlett  Packard  Model  200  C is  used  to  excite  the  primary  coils  of 
the  differential  transformers.  The  output  signal  frequency  is  set  at  7300  cps, 
this  value  being  determined  by  experiments  with  the  transformer  to  afford  good 
stability  and  sensitivity. 

2.  Power  Amplifier! 

It  was  found  necessary  to  supplement  the  oscillator  with  a power  ampli- 
fier due  to  the  low  impedance  of  6 differential  gages  in  parallel.  The  input  to 
the  power  amplifier  was  adjusted  so  that  2.3  volts  were  present  at  its  output 
terminals  when  fully  loaded.  See  Figure  19  for  circuit  details. 

3.  Amplifiers! 

A separate  amplifier  serves  each  pressure  cell.  The  electrical  signal 
taken  from  the  output  terminals  of  the  linear  differential  transformer  is  ampli- 
fied and  detected  by  means  of  the  circuit  indicated  in  Figure  19.  Vi  and  V2 
represent  straight  amplification  of  the  input  signal.  The  cathode  resistor*  of  Vi 
was  left  unbypassed  to  introduce  a small  amount  of  negative  feedback  for  stability 
purposes.  C,  was  found  necessary  to  reduce  a small  amount  of  parasitic  osillation 
that  occurred  from  feedback  from  associated  channels  wired  on  the  same  chassis.  V3 
acts  as  a modified  infinite  impedance  detector)  C„  is  the  detection  capacitor  act- 
ing in  conjunction  with  R^  and  R^^.  The  left  half  of  V3  is  biased  to  cutoff  by 
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proper  selection  of  the  grid  to  ground  voltage  of  the  right  half  of  this  tube. 
is  a difference  amplifier  providing  a high  current  source  for  the  low  impedence 
of  the  oscillograph.  Si  and  Ji  permit  the  introduction  of  a calibration  signal  to 
the  oscillograph,  which  is  normally  Inserted  between  the  cathodes  of  by  means 
of  the  jack  Ja . r^o  permits  the  selection  of  the  d.c.  level  of  the  output  signal, 
while  R17  permits  a fine  control  of  the  output  d.c.  signal  amplitude.  The  ampli- 
fier is  designed  for  a 7500  cps  carrier  signal  end  the  output  to  * 5/1  load  is 
linear  between  -8  to  +8  ma.  Two  such  amplifiers  were  constructed  on  a single 
chassis. 

4.  Power  Supplies! 

Two  300  volt,  70  ma  electronically  regulated  power  supplies  were  used 
in  conjunction  with  the  above  mentioned  amplifiers.  During  the  experimental  work 
described  in  this  report,  one  power  supply  energized  4 amplifiers,  while  the  other 
supplied  power  for  the  remaining  two.  The  Power  Supply  circuit  is  indicated  by 
Figure  19. 

5.  Recording  Oscillograph! 

The  Hathaway  Type  S-8B  recording  oscillograph  is  used  to  photograph  the 
output  traces.  A light  beam  is  focused  on  a mirror  set  on  the  mnvable  suspension 
of  a galvanometer  (No.  OA-2)  having  a natural  frequency  of  750  cps.  The  input 
signal  sets  the  galvanometer  suspension  in  motion,  causing  the  reflected  light 
beam  to  vary  accordingly  over  the  surface  of  a moving  photosensitized  paper.  The 
time  scale  of  the  record  is  obtained  by  a tuning  fork  arrangement  that  operates  an 
independent  shutter  and  light  source  system.  This  produces  a set  of  equally  spaced 
lines  on  the  photographic  paper.  The  overall  system  works  in  the  following  manner! 

The  oscillator  and  power  amplifier  energize  the  differential  transformer 
with  a carrier  frequency  of  7500  cps.  The  differential  output  of  the  secondary 
of  this  transformer  varies  according  to  its  core  position,  which  is  determined  by 
the  instantaneous  pressure  in  the  system  to  be  measurin'.  This  amplitude  modulated 
signal  is  then  amplified  and  detected.  The  detected  voltage  is  placed  upon  the 
terminals  of  the  oscillograph  galvanometers,  which  in  turn  determine  the  trace  of 
a light  beam  on  the  photo  sensitized  paper.  The  paper  is  developed  with  ordinary 
darkroom  techniques  and  sample  records  are  shown  in  Figure  8. 

H.  PgrfQrmaTKg 

Desirable  characteristics  of  a pressure  cell  are  reliability  of  indica- 
tion, ruggedness,  and  durability. 

1.  Reliability  of  Indication: 

In  general,  reliability  of  indication  means  that  once  a gage  is  calibrated, 
a given  indication  or  signal  represents  the  same  pressure  which  gave  that  indica- 
tion during  calibration.  Factors  tending  to  reduce  reliability  of  an  electronic 
diaphragn  gage  ares 

a.  Hysteresis 

b»  Shift  of  zero  position 
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Amplifier  for  Pressure  Cells 


Oscillator  Power  Amplifier 


toon 


Fig.  19:  Circuit  Diagrams  of  Electrical  Components  for 

Differential  Pressure  Cells 
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c.  Change  in  spring  constant  of  diaphragm 

d.  Low  frequency  response  (inability  of  diaphragm  and  sensing  element 
to  follow  varying  pressures) 

e.  Instability  of  the  electronic  element  under  conditions  of  normal 

use. 

Hvstareslsi  If  one  carefully  plots  the  output  curve  of  a diaphragm 
pressure  cell  Us  he  gradually  increases  the  pressure  and  again  as  he  gradually 
decreases  the  pressure,  he  usually  finds  two  curves  which  may  close  at  both  ends 
but  do  not  tough  in  the  central  region.  This  phenomenon  is  known  as  hysteresis 
and  is  due  to  failure  of  the  diaphragm  to  recover  immediately.  If  sufficient  tisie 
is  allowed  between  readings  for  the  diaphragm  to  recover  completely,  this  effect 
may  go  undetected.  Hysteresis  can  be  reduced  by  prudent  choice  of  diaphragm 
material  and  proper  mounting.  The  corrugated  beryllium  copper  and  Inconel  dia- 
phragas  used  in  the  gages  under  discucsicr.  showed  hysteresis  to  be  less  than  l£ 
of  maximum  applied  pressure*. 

Shift  in  Zero  Position*  This  effect  is  usually  due  to  insufficient 
clamping  at  the  edges  or  to  over  stressing  the  diaphragm*  In  the  case  of  flat 
diaphragna,  a considerable  amount  of  radial  tension  is  induced  itfien  the  diaphragm 
ia  deflected.  It  has  often  been  found,  especially  for  stiff  diaphragms,  that 
the  diaphragm  slips  away  from  its  clamping  at  the  edges  under  loading  and  then 
fails  to  slip  back  when  the  load  is  removed.  In  order  to  avoid  this  problem, 
diaphragms  have  sometimes  been  machined  out  of  a heavy  piece  of  metal,  leaving  a 
substantial  ring  at  the  edges  for  mounting.  One  successful  method  developed  on 
this  project  for  another  gage  (see  Reference  2)  Consists  of  sweating  a beveled  ring 
to  the  diaphragm  and  clamping  with  a conical  mating  ring  to  give  radial  tension. 

In  the  differential  pressure  cells  described  herein,  corrugated  diaphragms  soldered 
to  a heavy  ring  with  initial  radial  tension  show  no  zero  shift. 

Change  in  Spring  Constant  of  Diaphragm  This  effect  is  usually  due  to 
temperature  changes  and  for  all  practical  purposes  is  negligible.  When  used  In 
warm  or  cold  fluids,  the  gages  should  be  calibrated  at  the  same  temperatures. 

Frequency  Response!  If  a simple  spring  mass  system  is  excited  by  a 
periodic  <02  * of  the  same  frequency  as  its  natural  frequency,  the  force  is  said 
to  be  in  resonance  with  the  system,  and  the  latter  will  vibrate  with  a large; 
amplitude  than  the  exciting  force.  If  the  exciting  frequency  is  very  high  compared 
to  the  natural  frequency,  the  system  will  vibrate  at  its  natural  frequency.  If, 
however,  the  exciting  frequency  is  very  low  compared  to  the  natural  frequency,  the 
system  rill  vibrate  at  the  frequency  of  the  exciting  force.  In  obtaining  dynamic 
pressure  measurements  where  the  diaphragm  is  the  spring  and  the  mass  includes 
that  of  the  diaphragm  and  core  (see  Figure  17)  as  well  as  the  water  In  the  lead 
lines,  only  the  third  case  can  possibly  give  true  indications  of  pressure  pul3es. 

It  is  customary  to  assume  a gage  suitable  for  measuring  pulses  of  a frequency  up 
to  approximately  1/3  the  natural  frequency  of  the  diaphragm  sensing  element. 

Common  practice  Is  to  rate  pressure  cells  in  air  without  accounting  for 
the  mass  of  fluid  to  be  moved  in  lead  lines,  which  is  a variable  of  application. 

The  natural  frequencies  of  the  four  types  of  diaphragms  used  in  this  work  were 
determined  by  exciting  each  with  a 15  watt  speaker  driven  by  an  oscillator  and 
audio  amplifier.  The  normal  output  was  noted  on  an  oscilloscope  screen  as  shown  in 


Figures  20,  21.  nru»r.  resonant  frequency  was  applied,  tho  diaphragm  vibrated  with 
large  amplitude  and  the  signal  appeared  a&  in  Figure  20  b.  The  results  are  tabula- 
ted below* 


PUahMfB  AlMBhlY  Natural  Frequency  in  Air  fiflUXt  BXMtJM 


Full  Sire  Be-cu 
Half  Sire  Be-cu 
Full  Size  Inconel 
Half  Size  Inconel 


ZOO  cps 
1100  cps 
3000  cps 
6200  cps 


165+  cps 
220.  cps 
300  cps 


Of  greater  interest  in  the  study  at  hand  than  the  natural  frequency  of 
the  gage  in  air  is  the  natural  frequency  of  the  gage  as  installed  in  the  systqau 
To  approximately  obtain  this  frequency,  the  gages  were  set  in  the  running  position 
and  deliberately  subjected  to  water  hammer  shocks  by  opening  and  closing  the  foot 
valve  (see  Section  III-C) . These  shocks  caused  the  diaphragms  and  water  in  the 
lead  lines  to  vibrate  and  results  were  recorded  on  the  oscillograph  as  shown  in 
Figure  21.  The  large  oscillations  are  due  to  the  natural  frequency  of  the  water- 
air  system,  and  the  small  oscillations  are  due  to  water  hammer  of  the  system  with 
the  diaphragm  acting  to  effect  a more  elastic  system.,  These  frequencies  repre- 
sent those  which  the  diaphragms  can  pick  up  and  are  below  the  actual  natural  fre- 
quencies. Inasmuch  as  the  full  size  Ba-cu  ulaphragps  are  used  only  for  level 
measurement  a and  axe  quite  fragile,  they  were  not  subjected  to  the  water  hammer 
treatment.  The  recorded  frequencies  of  the  entire  system  for  the  finer  three 
diaphragm9  arc  shown  in  the  right  hand  column  In  the  above  tabulation. 

Instability  of  Electronic  EiamAntai  The  transformer  elements  employed 
in  these  pressure  cells  appear  to  be  very  stable  over  periods  of  at  least  several 
months.  It  has  been  found,  however,  that  it  is  difficult  to  obtain  oscillators 
and  amplifiers  which  are  stable.  For  this  reason,  it  is  advisable  to  make  fre- 
quent calibrations  during  test  runs. 

Ruooadnassi  The  gages  do  not  need  to  be  handled  with  extreme  care,  and 
only  i&en  the  most  sensitive  diaphragm9  axe  used  with  extremely  high  gain  settings 
are  the  cells  affected  by  nomal  room  vibrations* 

Durability*  Once  properly  waterproofed,  a cell  gives  long  trouble  free 
service  under  water.  In  seme  cases  they  have  not  been  opened  fox  adjustment  for 
several  months* 

f.  ARPUcaUon 

1.  Selection  of  Diaphrag/ii-Transrorraer  Combination* 

Selection  of  the  best  diaphragm-transformer  combinations  for  a given 
differential  pressure  measurement  is  largely  a matter  of  experience.  Of  course., 
for  pressures  over  50  psi,  only  the  "half  size  Inconel"  diaphragm  is  strong  enough 
to  resist  permanent  set,  while  the  beryllium  copper  diaphragms  should  not  be  used 
for  pressures  of  more  than  about  12  ft.  and  20  ft.  of  water  for  the  full  and  half 
sizes  respective’  y.  However,  even  the  full  size  Inconel  diaphragms  can  ba  used  »dth 
fair  success  in  ranges  of  0 to  10  ft.  water  for  full  10"  chart  deflections,  and 
the  half  size  Inconel,  good  for  about  100  psi,  has  been  employed  to  give  full  10" 
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chart  deflections  fez  zanies  as  lew  as  0 to  15  psi.  It  Is  thus  apparent  that 
there  is  considerable  overlapping  in  adaptability* 

It  has  been  observed  that  these  corrugated  diaphragns  can  be  deflected 
to  about  of  their  free  diameters  without  Incurring  any  residual  set.  For  the 
smallest  diaphragp,  the  half  size  Inconel,  which  is  0.5"  diameter,  this  would 
mean  a deflection  of  0.020",  or  the  full  range  of  the  largest  transformers  In  use. 
Hence,  It  would  not  be  possible  to  overload  any  diaphragm  without  noting  the  devia- 
tion from  linearity  of  the  output.  From  the  standpoint  of  obtaining  instantaneous 
measurements,  the  smallest  total  displacement  Is  the  most  desirable,  since  this 
requires  the  smallest  quantity  of  water  to  be  set  in  motion  in  the  lead  lines.  On 
the  other  hand,  a vory  small  movement  of  the  diaphragp  for  full  chart  deflection 
requires  the  electrical  amplification  of  the  signal  to  be  highj  this  often  results 
in  a shift  in  calibration  due  to  only  minute  changes  in  the  primary  input  voltage. 
Experience  has  shown  that  a diaphragm  transformer  combination  tdiich  under  full 
pressure  ranae  allows  the  diaphxagpi  to  travel  about  6CJ6  of  the  rated  range  of  the 
transformer  Ci.e.  0.003"  deflection  of  diaphragm  for  a 0.005;  transformer),  give® 
very  satisfactory  results.  It  should  be  noted  th«.t  the  full  chart  deflection  can 
be  obtained  for  ranges  other  than  those  starting  from  zero  (viz.,  3 psi  to  9 pal 
rather  than  0 to  6) • 

2.  Setting  Null  Position: 

Characteristics  of  the  differential  transformer  used  in  the  differential 
pressure  cell  are  such  that  the  output  is  linear  with  travel  of  the  core  within 
the  range  of  the  instrument.  When  the  core  is  centered  between  the  two  secondary 
coils.  It  is  at  a null,  or  zero  position.  The  net  output  of  the  transformer  here 
la  zero.  The  signal  increases  as  the  core  moves  away  from  this  null  point  iii  either 
direction;  however,  in  the  vicinity  of  the  null,  the  output  changes  only  s lightly 
with  core  movement.  For  this  reason,  the  transformer  is  usually  operated  in  a 
range  which  starts  with  the  core  a small  distance  from  the  null  and  allows  it  to 
proceed  in  the  same  direction  from  the  null  with  increasing  dlaphragp  deflection. 

Of  course,  comparable  results  can  be  obtained  by  starting  some  distance  from  the 
null  point  and  allowing  the  core  to  approach  the  null  position,  so  long  as  the 
"dead"  region  in  the  vicinity  of  the  null  is  avoided.  In  the  latter  case,  the  sig- 
nal would  decrease  with  increasing  deflections. 

The  procedure  for  setting  the  null  position  is  best  explained  by  referring 
to  Figure  17.  The  cover  is  removed  from  the  gage  and  the  assembly  of  diaphragm, 
cage,  and  adjusting  nut  are  removed  from  the  gage  body.  The  primary  coil  is 
excitSd  ( normally  7500  cps  and  about  2.5  volts),  and  the  output  signal  is  observed 
as  a sine  wave  on  an  oscilloscope  screen.  By  turning  the  adjusting  nut  in  or  out 
of  the  cage,  the  core  i3  moved  relative  to  the  transformer  coils.  The  null  point 
is  reached  when  the  signal  on  the  oscilloscope  screen  Is  a minimum.  In  practice, 
the  nut  is  turned  sufficiently  to  permit  the  core  to  advance  slightly  beyond  the 
null  point,  so  that  when  the  diaphracp  is  loaded  the  motion  will  continue  In  the 
same  direction. 

When  the  desired  setting  is  attained,  the  lock  nut  is  tightened,  the 
assembly  pushed  back  into  the  gage  body,  and  the  cover  replaced. 
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3.  Procedure  for  Filling  Gages  with  Water  and  Bleeding  Piezometer  Tape  and 
Lead  Linesi 

No  air  pockets  can  be  trapped  in  the  gage  or  lead  lines  if  the  recorded 
pressure  indications  are  to  be  "true".  Air  pockets  compress  with  increased  pres- 
sures and  cause  wore  flow  to  be  set  up  in  the  lead  lines  than  that  required  to 
depress  the  diaphraga.  The  inertia  of  the  water  thus  set  in  motion  induces  errors 
both  in  the  magnitude  of  the  pressure  changes  and  in  the  time  it  which  the  changes 
occur.  Air  pockets  are  reduced  by  bleeding  the  lead  lines  and  piezometer  taps  and 
by  evacuating  the  pressure  cells  before  filling. 

The  procedures  for  evacuating  the  pressure  gages  are  best  explained 
by  referring  to  Figure  9t 

a.  A Bourdon-type  vacuum  gage  is  attached  at  valve  #4. 

b.  A short  length  of  copper  tubing  is  attached  at  valve #2t 

c.  A chsmlcal  laboratory  type  aspirator  is  attadied  to  a faucet. 

Rubber  tubing  leads  from  the  vacuum  side  of  the  aspirator  to  the 
above  short  length  of  copper  tubing.  This  juncture  is  then  im- 
mersed in  a Jar  of  water. 

d.  All  valves  on  the  gage  "Christmas  tree"  axe  opened  except  #1. 

e.  When  the  vacuun  reaches  28.3”  o.  tg.  or  more*  the  rubber  tube  is 
moved  from  the  short  copper  tubing  and  water  enters  the  cavity. 

When  the  gage  is  filled,  it  is  placed  on  the  supporting  rack  and  the 
pressure  lead  lines  are  attached.  Each  is  then  bledrthrough  the  gage  valves  for 
several  minutes  before  each  set  of  runs.  A piece  of  tubing  attached  to  the  bleed 
valve  (No.  l)  facilitates  this  bleeding  by  directing  the  bled  water  into  the 
bottom  tank. 
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